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We investigate the generation of multipartite entangled state in a system of A^ quantum dots 
embedded in a microcavity and examine the emergence of genuine multipartite entanglement by 
three different characterizations of entanglement. At certain times of dynamical evolution one can 
generate multipartite entangled coherent exciton states or multiqubit W states by initially preparing 
the cavity field in a superposition of coherent states or the Fock state with one photon, respectively. 
Finally we study environmental effects on multipartite entanglement generation and find that the 
decay rate for the entanglement is proportional to the number of excitons. 
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I. INTRODUCTION 

Quantum information processing offers important ca- 
pabilities for quantum communications and quantum 
computation in a variety of physical systems, and the 
solid-state implementation is one of the most promis- 
ing candidates. Semiconductor quantum dots (QDs) are 
appealing for the realization of quantum computer and 
quantum information processing. Spatial confinement of 
excitons in three dimensions leads to discrete energy lev- 
els, and the main causes of phase decoherence, namely 
scattering events, are highly suppressed 0. In this re- 
spect, QDs are very promising compared to other semi- 
conductor structures with band states. Moreover a quan- 
tum computer scheme based on quantum dot arrays is 
scalable up to > 100 qubits. 

Entanglement is a essential source for quantum com- 
puting and quantum information processing. It is well- 
known that a controlled-NOT gate can generate a max- 
imally entangled state. Various schemes for the realiza- 
tion of controlled-NOT gates in QDs [^|j3|Q] are proposed 
with either the electron spin ^|JJ|] or the discrete elec- 
tronic charge degree of freedom as qubits For 
example, Imamoglu et al. Q proposed a scheme that re- 
alized a controlled-NOT gate C^ OT between two distant 
quantum dots i and j via the cavity quantum electro- 
magnetic dynamics (CQED) techniques. Here the first 
superscript i in C^ OT denotes the control and the second 
denotes the target. These controlled-NOT gates in turn 
can be used to generate a tripartite maximally entangled 
Greenberger-Horne-Zeilinger (GHZ) M state as |0)<g>|0)<g> 
|0) + |1) ® |1) ® |1) = a 3 O T^ 2 o T (|0) + |1» ® |0) ® |0>. 
Quiroga et al. showed [[l0| that an optically controlled ex- 
citon transfer process can generate maximally entangled 
Bell states [jllj and GHZ states. To demonstrate that 
excitons in coupled quantum dots are ideal candidates 
for reliable preparation of entangled states in solid-state 
systems, Chen et al. |l2| experimentally investigated how 
to optically induce entanglement of excitons in a single 
gallium arsenide QD: the entanglement is identified by 



the spectrum of the phase-sensitive homodyne-detected 
coherent nonlinear optical response. 

Although no experimental observation of entanglement 
of many-exciton systems has yet occurred, the theoreti- 
cal studies are necessary if we desire to have a large-scale 
quantum computing device. In a recent investigation by 
Liu et al for generating bipartite entangled coherent 
exciton states in a system of two coupled quantum dots 
and CQED with dilute excitons, they found that the bi- 
partite maximally entangled coherent exciton states can 
be generated when the initial cavity field is in an odd 
coherent state. The relation between entanglement of 
exciton and photon numbers in the cavity was also dis- 
cussed. 

It is interesting to extend to multipartite entanglement 
in such systems. Multipartite entanglement is not only 
of intrinsic interest itself but also of practical importance 
in quantum information proposals such as quantum tele- 
portation and quantum cryptography. One motivation to 
consider multipartite entanglement in many QDs is that 
around 10 QDs can be embedded in a microdisk struc- 
ture and coupled to a single cavity mode in the quantum 
information process scheme proposed by Imamoglu et al. 
We will first extend the model of Ref. [11] to the case of 
many QDs , and then study how to generate multipartite 
entangled states and examine multipartite entanglement 
in such systems. 

The paper is organized as follows. In Sec. II, we present 
a model for coupled QDs in CQED and determine the ex- 
act state vector at any time t, which will be shown to be 
a multipartite entangled state in some time ranges. In 
Sec. Ill, we use three different methods to examine the 
multipartite entanglement in the state. Then in Sec. IV, 
we show that multipartite entangled coherent exciton 
states can be generated at certain times during the dy- 
namical evolution. In Sec. V, we analyse the effects of 
environment on the generation of multipartite entangle- 
ment. The conclusions are given in Sec. VI. 
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II. MULTIPARTITE ENTANGLED STATES 

We consider quantum dots embedded in a single-mode 
cavity. We assume that the QDs are sufficiently large and 
there are only a few electrons excited from valence band 
to conduction band |13|. In the assumption of low exci- 
ton excitation density, exciton operators can be approx- 
imated by boson operators, and all nonlinear dynamics 
such as exciton-exciton interaction can be neglected. As 
the distance between any two excitons is also assumed to 
be large, the interaction between any two excitons can 
be safely neglected. There exists a resonant interaction 
between the excitons and the cavity field. The Hamilto- 
nian under the rotating wave approximation is given by 

(h = l) 



H 



cubtbo 



N 



N 



Y,g n (btbo + bib n ), (1) 



where b (bo) is the creation (annihilation) operator of 
the cavity field with frequency ui, b^bn) [n = 1, 2, • ■ • , N) 
denote the creation (annihilation) operator of the n th ex- 
citon with the same frequency u>, g n is the coupling con- 
stant between the cavity field and n th exciton, and N is 
the number of excitons. 

It is convenient to write the coupling constants gi in 
generalized spherical coordinates Jl5| , ^6[ as 

5i = G cos 0i, 

.92 = G sin tfii cos 02 > 

93 = G sin </>i sin 02 cos 03 , 

gN-i = G sin 4>i sin 02 • ■ ■ cos0/v_i, 
g n = G sin 0i sin 02 ■ ■ ■ sin0jv-i, 



where G = yX^fcLi 9 k 2 - The evolution operator corre- 
sponding to the Hamiltonian is then obtained as 

U(t) = ' ' ' ^3(02)^2(01)^0,1(0^1,2(01) 

x\/ 2)3 (0 2 ) • • • V N -x, N {^ N -{)U Q {t), (2) 



where 



7 iJ (0 i )=ex P [0 i (&t6 i -6t 6i )] (i^fej), 
V 0A (t) =cxp[-iGt(blb 1 + b\b )}, 



N 



U (t) = exp[-iwt(6j6o + Y] b\h)} 



=1 



(3) 



Both Vtj((j>i) and Vb,i(t) are operators for the beam split- 
ters. We see that the evolution operator can be written 
as a free evolution Uo(t) followed by a series of beam 
splitter operators. 

Let us assume that the initial state is \ip(0)) — \a) <8> 
|0) 1 ® • • -(g) \0) N , i.e., the cavity field is in a coherent state 



\a) a with a^O and all the excitons are in vacuum states. 
After the action of the unitary operator U(t), the state 
evolves into 



|V>(*)) = \<*o{t), ai(t), • • • , a N (t)) ...N, 

where \a (t), • • • , ajy(*))o-iv = \<*o{t))o ® |ai(i))i ( 
\a N (t)) N and 



(4) 



0*0 (*) = acos{Gt)e~ tu> \ 
a n (t) = -iasm(Gt)g n /Ge- tuJt 



,n=l,2,.-.,JV. 



(5) 



Thus, if the initial state of the cavity field is in a super- 
position of coherent states, the resulting state will be an 
entangled coherent state. Let us assume the initial state 
of the cavity field to be in a superposition of two coherent 
states | ±a) and the excitons to be in vacuum states, i.e., 

|*(0)) = [2 + 2cos6>exp(-2H 2 )n 1/2 

(Mo+e^l-aUslOhB-.-SlO)*. (6) 

Specifically, for 9 = 0, n, and 7r/2, the cavity superpo- 
sition state reduces to even, odd fl7|] , and Yurke-Stoler 
]]~8|] coherent states, respectively. Then the state vector 
at time t is given by 

|*(t)) - 7V(M*), ai(i), • • • , ajv(t))o...iV 

+e ie | - ao(t), -ai(t), • • • , -ajv(t))o-iv), (7) 



where 



JV 



AA= [2 + 2cos6']^[p fe (<)]- 1 



-1/2 



(8) 



fe=0 



is the 
normalization constant, and p k (t) = exp(— 2 1 cka; (i) | 2 ) is 
the overlap of the two coherent states | ± a k (t)). The re- 
sulting state is a multipartite entangled coherent 
state 

We now choose the orthogonal basis [|o[ 
\0) k = \a k (t)), \l) k = (\-a k )-p k (t)\0))/M k (t), (9) 



where M k (t) = y^— Pk(t) 2 - It then follows that 

\-a k (t)) =M k (t)\l)+ Pk (t)\0). (10) 

Using this basis the state vector at time t can be rewritten 
as 

|tf (t)) = AA{|0) ® • •• ® |0)jv + e ld [Mo(t)\l) + po(*)l°)o] 
<g> •• • ® + (t) 1 0) at]}. (11) 

After the 'encoding' this state is a multiqubit state. Then 
we can fully exploit the sophisticated tools available for 
examining multipartite entanglement of qubits to study 
our state and determine if it is genuinely multipartite 
entangled. 
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III. EXAMINATION OF MULTIPARTITE 
ENTANGLEMENT 



tion fidelity |23j 
concurrence 




We study the multipartite entanglement by examining 
(i) Mermin-Klyshko inequality |2l],|22j , (ii) state prepara- 
and (iii) the square of the multiqubit 
The first two, namely the Mermin- 
Klyshko inequality and the state preparation fidelity, are 
related to two sufficient conditions that distinguish be- 
tween genuinely TV-partite entangled states and those in 
which only M particle are entangled (M < TV). The 
last one, the square of the multiqubit concurrence, is not 
by itself a measure of TV-particle entanglement, but it 
appears to be related to some kind of multipartite entan- 
glement. That is, we can gain some information about 
the degree of multipartite entanglement by calculating 
the square of the multiqubit concurrence. 



A. Mermin-Klyshko inequality 



B(p) = 



\(B N )\-2 N / 2 

2(W+l)/2 _ 2 A72' 



(17) 



Then the state p is TV-partite entangled when B(p) > 
and maximally entangled when B(p) = 1. 

Let Ai — a x and A\ = a y for any i. Then the Bell 
operator Bn and B' N become Q 

B N = 2( N + 1 )/ 2 ( e -^~ ( 7f v + e^a®"), 

B' N = ^ N+1 ^ 2 (- le -^af N + ie^a® N ), (18) 

where f3 N = tt/4(TV - l),a+ = |0)(1|, and cr_ = |1)(0|. 
For our state \^(t)} (|ll|) we have 



B(|*(*))<*(t)|) 

2 (iV+l)/2 



1 + cos6»e- 2 l Q l 2 

N 

X 

k=0 



cos(6 — j3 N ) + cos(/37v)e' 



-2\a\' 



nv 1 



pi- 



rn 



Let us first use the Mermin-Klyshko inequality [£1] 22| 
to examine the TV-partite entanglement. This inequality 
generalizes the Bell inequality |TTJ and Clauser-Horne- 
Shimony-Holt (CHSH) inequality j27j, which not only 
tests the predictions of quantum mechanics against those 
of local hidden variable theory but also distinguishes en- 
tangled from non-entangled states [ p8|p4| ] . The Mermin- 
Klyshko inequality is pip^ ] 

l(MI<2, (12) 
where Bn is the Bell operator defined recursively as 

B N = ^{A n + A' n ) <gi B N -i + ^{A n - A' n ) <8> B' N _ X , 

(13) 

B' N is obtained from 23 at by exchanging primed and un- 
primed terms, 

B'n = \{A N + A'n) ® B' N ^ -\{An- A' n ) ® Ejv-i, 

(14) 

Bi = 2Ai, and B[ = 2A[. All A t and A[ are dichotomous 
observables. Let S n denote the set of all TV-particle states 
and S^~ x denote the subset of those states which are at 
most (TV — l)-partite entangled. Then, from the results 
of Refs. [p8|j24|], for a state p we have 



Vp eSj^- 1 : |<0„>| <2*l\ 
VpeS N : \(B N )\ < 2 ( N+1 )/ 2 , 



(15) 
(16) 



which implies that a sufficient condition for TV-partite 
entanglement is the violation of the inequality given by 
Eq. dlq). Now we define a quantity B 



We choose 9 = (3n- Then the above equation reduces to 

N 



B=2 (^)/ 2 nv 1 -^ 

k=0 



Therefore this sufficient condition becomes 



A' 



l[Jl-pl>l/V2. 



(20) 



(21) 



Numerical results are provided later. Next we discuss 
another method based on state preparation fidelity to 
examine multipartite entanglement. 



B. State preparation fidelity 

The so-called state preparation fidelity T of an N- 
qubit state p is defined as 



F(p) = (V'GhzIpIV'Ghz), 



where 



IV^GHZ) = ^=(|0) 



10) + e^ll) 



(22) 



|1» (23) 



is the GHZ state. A sufficient condition for TV-partite 
entanglement is given by P3,E4| 



Hp) > i/2. 



(24) 



From Eqs. (|22|) and (|llj) the state preparation fidelity 
for the multipartite state \^(t)) is obtained as 
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*W)X*(*)I) 
1 



N 



4 + 4cos 6e- 2 \ a \ 2 



i+n( i ~^)+ e_4H 



fe=0 



iV 



+2cos6»e- 2|Q| +2cos(6»-7) JJ Jl 



Pt 



k=0 



N 



k=0 



(25) 



Let 9 = 7 = 7r/2. The above equation then reduces to 



(jv \ • 

fc=0 / 



/4 + e- 4 l Q l/4, (26) 



and the sufficient condition ( p4| ) becomes 

/ AT ^ 2 



i+n v 1 -^ + e - 4H >2. 



fc=0 



For convenience wc define 



F = IT - 1. 



(27) 



(28) 



When F > the state is (N + l)-partite entangled. 

In Fig. 1 we plot the quantity B and F against time t. 
Henceforth we assume g% = gi = ... = and G = 1. 
The period of B and F with respect to time is 2tt. We can 
clearly see the two sufficient conditions can be satisfied in 
some time ranges which means that the state is a genuine 
multipartite entangled state. Moreover, we find the time 
range in which the sufficient condition is satisfied based 
on the state preparation fidelity is larger than that based 
on the Mermin-Klyshko inequality. 



C. Square of the multiqubit concurrence 

Recently Coffman et al |g5| used concurrence Q to 
examine three-qubit systems, and quantified the amount 
of tripartite entanglement in three-qubit systems by the 
quantity tq,i,2 25 



To, 1,2 — Cq(12) ~ Q)i ~ Q)2> 



(29) 



where Co(i2) denotes the concurrence between qubit 1 
and qubits 2 and 3. Applying the general result for con- 
currence of bipartite nonorthogonal pure states to 
the state (O) for N = 2 yields 



0(12) 



V(l-p 2 )(l-p 2 P 2 ) 

1 + P0P1P2 cos 9 



(30) 



By the standard method for calculating the concurrence 
we find 



01 



C02 — 



P2X /(i-p 2 )(i-p 2 ) 

1 + cos 6p a pip 2 
Plv /(1-P 2 )(l-^ 
1 + cos 9p a pip2 



(31) 



The concurrence C02 is obtained by a transformation 
1 < — ► 2 in the expression of Cq\. From Eqs. (p9|), (|3T 
and (|3l|), we obtain 



70,1,2 = 



n„= (i-^) 



(i 



cos ( 



; l«l 2 ) 2 



(32) 



Now consider the case of odd TV excitons (N > 2). 
Wong and Christensen J2(| proposed the square of the 
multiqubit concurrence as a potential measure of the mul- 
tipartite entanglement for an even-number pure qubit 
state. The (N + lVqubit concurrence for a pure state 
\ib) is defined as lEQl 



where a v 



-i(|0)(l| 



(33) 



Jl)(0|) is a Pauli matrix. 
Then applying Eq. ( |33| ) to the state (|ll|) leads to the 
square of the multiqubit concurrence 



nto(i-fn) 



r 0,l, 



C! 



0,1, 



,N 



(1 



cos f e 



-2|a| : 



) 2 



(34) 



From the expression of the square of multiqubit concur- 
rence we can see that the multiqubit concurrence is un- 
changed by permutation of qubits, which implies that it 
really represents (N + l)-partite entanglement. We note 
that although the above formula is obtained for odd N, by 
comparing ( |34| ) and (|32|) , it is also applicable to N = 2. 
We also note that Eqs. (|34j) and are applicable to the 
more general state (|^o) ® |"0i) ® • ® \^Pn) + e l0 \fio) <8> 
1 0i) <£>•■■<£> I 0jv)) with real overlap (V't|^t)(i = 0, 1, ...,N) 
up to a normalization constant. 

In Fig. 2 we plot the square of the multiqubit concur- 
rence against time for N = 2, N = 3, and N — 5. We 
observe that the entanglement periodically reaches its 
maximum twice per period, and the multipartite entan- 
glement is suppressed with the increase of the number 
of excitons. One way to overcome this suppression is to 
increase the parameter \a\ 2 . From Eq. (34) we know that 
the larger the parameter |a| 2 , the larger the multipartite 
entanglement. 

From Eq. (|34] ) we know that the square of the (N + 1)- 
qubit concurrence reaches a maximum value at 9 = it 
when other parameters are fixed. This implies that the 
best input state is the odd coherent state in order to 
generate multipartite entanglement. On the other hand 
we know that the average photon number of the cavity 
field in the initial state |^(0)) is 



(blb ) 



1 - cos6» e - 2 l a l 2 
1 + cos 6>e- 2 H 2 



(35) 
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We see that the mean photon number of the cavity field 
reaches its maximum when the field is in the odd coher- 
ent state (9 = 7r) when other parameters are fixed. The 
mean photon number represents the energy of the system. 
Therefore, it turns out that the more energy contained 
in the initial cavity field, the larger the multipartite en- 
tanglement. 

By three different methods we have examined multi- 
partite entanglement and found that the state ( pi] ) is a 
genuine multipartite entangled state over a large range 
of parameters by each of the indicators. 



IV. MULTIPARTITE ENTANGLED COHERENT 
EXCITON STATES 

It is clear that the cavity field and excitons are de- 
coupled when t — nir /(2g) for integer n. Without loss of 
generality we examine the state vector at times 7r 

and 7r/2 in one period. From the discussion of the last 
section, at these times there is no global multipartite en- 
tanglement. At time t = 7r, the state vector returns to 
the initial state (||) with a — > —ae~' luJV , and the exci- 
tons are in the vacuum state. It is interesting to see that 
at time t = tt/2 the cavity field is decoupled from the 
excitons; however the excitons are left in a multipartite 
entangled coherent exciton states given by 



I*) 



+e *| -ft-/?,.- 



■N 



-/3)i...n), 



(36) 



where (3 — —i(a/yN)e lW7r / 2 , and the normalization 
constant Af' = [2 + 2cos6>exp(-2|/3| 2 )]- 1/2 . Now the 
cavity field is in a vacuum state, and all its energy has 
been transferred to the excitons. 

Multipartite entanglement for the entangled coherent 
exciton state | ^excitons can be studied using the methods 
of the last section. From Eqs. @, @, and @ we ob- 
tain the quantities B, and the square of the multiqubit 
concurrence for our state I ^excitons as 



B = 2 N ' 2 {1- 



-A\a\ 2 /N^N/2 



(37) 



n, 



T = [1 + (1 - e -4|«| 2 /A r )A r /2]2 /4 + e -4M 2 /4j ( 38) 
H _ e -4|a| 2 /iV)iV 



(1 



cos tie 



-2\a\- 



(39) 



In Fig. 3 we plot the quantity B, F, and r against 
\a\. We observe that the state is multipartite entangled 
when \a\ is large enough. Quantitively the corresponding 
Bell inequality is violated when \a\ > 1. 601, and the 
state preparation fidelity T is larger than 1/2 when \a\ > 
1.228. We also observe that the square of the (N + 1)- 
qubit concurrence is significantly larger than zero only 
if | a | is large enough. For fixed N and very large \a\ 
the square of the (N + l)-qubit concurrence Ti,...,^ ~ 1 
which implies that the entangled coherent exciton states 
becomes a GHZ-likc state. On the other hand t\ ... n « 



for fixed N and small enough \a\ 2 . As discussed in 
Ref. |2(J the state |^ , ) cxc j tons with 8 = ir reduces to the 
multiqubit W state |3^] in the limit of \a\ 2 — > 0. It means 
that we can also prepare the W state in our system at 
time t = 7r/2 with the initial cavity field in a Fock state 
with one photon and all excitons initially in the vacuum 
states. For the case of only two excitons the W state 
is just the maximally entangled state (one Bell state) as 
discussed by Liu et al. fll3||. 



V. EFFECTS OF ENVIRONMENT ON 
MULTIPARTITE ENTANGLEMENT 

Environmental losses and decoherence are important 
effects in quantum information processing 33 1 . The life- 
time of both cavity photon and exciton is generally con- 
sidered to be on the order of picoseconds. Whereas, if 
we assume fig = 0.5 meV [ [l3"| , the time we need to get 
maximal entanglement in our model is also on the or- 
der of picoseconds. However, dynamical evolution suffers 
from decay of photons and excitons. Recent experiments 
show that the lifetime of photons and excitons can be pro- 
longed Q,|||j3^]. In particular, for cavity decay, recent 
experiments display an elongated decay time of photons 
in the microwave domain ]3q ]. Thus, we only consider 
the decay of excitons in the following discussions. 

We follow the method of Ref. [jl3| , and assume that the 
environment is at zero temperature and the system dis- 
sipates by interaction between excitons and a multimode 
electromagnetic field. Under the rotating-wave approxi- 
mation we write the Hamitonian as 



N 



g^2(blb + blb n ) 



71=1 

N 



^2^2 X k (bia k + b n c 



kJ' 



(40) 



n—l k 



where at (a k ) denotes the creation (annihilation) opera- 
tor of the multimode magnetic field with frequency uj k - 
We assume g\ = gi = ■ ■ ■ = g implying that the ex- 
citons are equally coupled to the cavity mode, and n- 
independence of the X k implying that the excitons are 
also equally coupled to the environment. 

Now we use the Heisenberg picture to study the prob- 
lem. The Heisenberg equations for related operators are 
obtained as follows: 

JV 

db Q /dt = iujb + ig^ bn ' 

n=l 

db n /dt = iub n + igb + i^X k a kl 
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N 



N 



da k /dt = iu k a k + i\ k ^ b n . 



(41) 



n=l 



Note that we use a slightly different Heisenberg pic- 
ture for our purpose to obtain the final state vec- 
tor at time t. We let an operator evolve as A(t) = 
cxp(—iHt)A(Q)cx-p(iHt) and it satisfies idA(t)/dt = 
[H,A(t)]. ^ 

Now we introduce the collective exciton operator b c = 
l/^Eti''" which can be considered as a boson since 
[b c , bl] = 1. In terms of the collective boson operator we 
rewrite the above equation as 

dbo/dt = iubo + igNb c , 

db c /dt = iub c + ignbo + i ^ \N,ka k , 

k 

da k /dt = iuj k a k + i\N,kb c , (42) 
where g^ — g^fN and \N,k = AfcV^V- If we let 

b Q = Boe**, b c = B c e^\ a k = A k e iut , (43) 
Eq. (H|) reduces to 

dBo/dt = ig N B c , 



dB c /dt = ig N B + i ^ \N,kA k , 

k 

dA k /dt = i(uj k - u)A k + i\ N , k B c . 



(44) 



To solve the above equation we make the Laplace trans- 
form 



f(s) = / e- st f(t)dt. 
Jo 

We obtain after minor algebra 
sB =B Q {0)+ig N B c , 



(45) 



(46) 



X N ,kA k (0) 



BJO) + ig N B + i V — 

B c = —— — ^ _ (4?) 

S + N J2k s+i(J° k -u) 

The above equation cannot be solved exactly. So we re- 
sort to the Wigner-Weisskopff approximation (37). After 
the approximation, from the above two equations we ob- 
tain 



Bn 



(s + ^)Bq(O) + l9 VNB c (0) - E fc a :+^^ 



NT , 



Ngi 



(48) 



where T = 27re(w)A 2 (u;) and e(ui) is a distribution func- 
tion of the multimode electromagnetic field. We assumed 
that Auj — — J duj k e(u) k )\ 2 (us k )/(u) k — u>) — in the 
derivation of the above equation. 

From Eq. ( f48| ) we obtain the operator bo{t) in the 
Heisenberg representation 



b (t) = u(t)6o(0) + V(t) b n(0) + «»fc(*)Ofc(0) s 



71=1 



u(t) = e- Nrt/4 [cos(A N t) + —— sin(A Jv t)]e i ' 
v{t) = - — e ' sm{A N t)e , 



(49) 



where A w = \/Ng 2 - iV 2 r 2 /16. Let the cavity field be 
an odd coherent state and other systems in the vacuum 
states. Then, from Eq. ( f49| ) we obtain the state vector at 
time t as 

|*(t)) = [2-2exp(-2M 2 )] 

x(|a u )® \a v )® N ® \[\a w , k ) 

k 

-\-a u )®\-a v )® N ®\{\-a w ^)), (50) 

k 

where a u — au*(t), a v — av*(t), and a Wik — aw k {t). 

We use the state preparation fidelity J- to examine 
multipartite entanglement in the above state and choose 
the following GHZ state for consideration, 

I^ghz) = ^(K> ® \a v )® N - \ai) ® \air N ), (51) 



where 



\a£) = (I - a x ) - p x \a x ))/M 3 
p x = exp(-2|a x | 2 ), 

M x = \Jl- pi, x = u,v. 



(52) 



From Eqs. @) and (|T|) we obtain the state preparation 
fidelity as 

^=(*(t)|^GHZ><VGHz|*(t)) 
1 



4 _ 4e -2|a| 
+ (M U M? -puP^ f 



[1 + 2e- 2 l Q l {M u M^p- l p- N - 1) 

(53) 



In the derivation of the above equation we have used the 
relation 



H 2 = \a u \ 2 + N\a v \ 2 +Y,\ a 



w.k | 



(54) 



which results from the normalization of the state ^(t)) 
and implies the energy conservation. 

Figure 4 shows the fidelity against time t for different 
numbers of excitons. From the figure we see that the fi- 
delity is not a periodic function of t due to the dissipation 
of energy to the environment. When there is no dissipa- 
tion (r = 0), the fidelity attains the ideal case of being a 
periodic function of t. We observe that multipartite en- 
tanglement occurs only in the beginning of the evolution. 
When the number of excitons becomes larger, the gen- 
eration of the multipartite entanglement becomes more 



G 



difficult. In the limit of t — > oo the fidelity becomes 0.5 as 
we expected. In this case there is no multipartite entan- 
glement and all the energy of the cavity-excitons system 
dissipate to the environment. The environment dimin- 
ishes the generation of multipartite entanglement when 
the number of excitons increases. From Eq. ( [i"9|) we can 
see that the decay rate is proportional to the number of 
excitons when the intensity of the cavity field is fixed. In 
addition, we find a similar result as that of Ref. [jl3| if 
we fix r and vary \a\, i.e., the multipartite entanglement 
decays rapidly with increasing the cavity field density. 



VI. CONCLUSIONS 

In conclusion, we have investigated the dynamical evo- 
lution of multipartite entanglement in a system of quan- 
tum dots embedded in a microcavity. The entanglement 
is studied via two sufficient conditions for multipartite 
entanglement and the square of the multiqubit concur- 
rence. We observed the global multipartite entanglement 
and at certain times the entanglement becomes maximal. 
We can also produce the multipartite entangled coher- 
ent exciton states and multiqubit W state by preparing 
different initial states. Finally, we study the effects of 
environment on the generation of multipartite entangled 
states, and find that the decay rate is proportional to the 
number of excitons. We also find that the entanglement 
decays rapidly with increasing the cavity field density. 

Although multipartite entanglement studied here has 
not yet been observed experimentally, the potential ap- 
plication of excitons in quantum computing as well as 
rapid development of CQED technique suggest that our 
analysis will find applications in this field. 
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FIG. 1. The quantities B (solid line) and F (dashed line) 
jainst time. For B we choose 6 — /3n- For F we choose 
= 7 = 7r/2. The parameter G = 1, |a| 2 = 3, and N = 3. 




FIG. 2. The square of the multiqubit concurrence against 
time. The parameter g = 1, |aj 2 = 0.9, 9 — it, and all gi 
are equal. The solid, dashed and dotted lines correspond to 
N — 2, 3 and 5 respectively. 
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FIG. 3. The quantity B (solid line), F (dashed line), and 
square of the multiqubit concurrence (dotted line) against \a\. 
The parameter N — 5. 
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FIG. 4. The quantity F — 2T — 1 for different numbers of 
excitons: N = 2 (solid line), N = 3 (dashed line) and iV = 4 
(dotted line). The parameters \a\ 2 =3,g =1, and T = 0.5. 
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